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I discuss the structure and composition of massive (two solar-mass) neu- 
tron stars containing hypernuclear and deconfined quark matter in color 
superconducting states. Stable configurations featuring such matter are 
obtained if the equation of state of hadronic matter is stiff above the satu- 
ration density, the transition to quark matter takes place at a few times the 
nuclear saturation density, and the repulsive vector interactions in quark 
matter are substantial. I also discuss our recent progress in understand- 
ing the cooling of massive compact stars with color superconducting quark 
cores. 

PACS numbers: 97.60. Jd, 26.60.Kp. 95.30.Sf 

1. Introduction 

The heavy-ion collision experiments and compact stars provide mutu- 
ally complementary channels to address some of the outstanding challenges 
of the modern particle and nuclear physics. While the energy densities 
achieved in these experiment are overlapping with the range that must ex- 
ist in neutron (or more generally compact) stars, their astrophysical studies 
still provide a window on the properties of dense matter that may be difficult 
or even impossible to extract otherwise. 

The masses of neutron stars are the most sensitive among their inte- 
gral parameters to the equation of state (hereafter EOS) at high densi- 
ties. Therefore, pulsar mass measurements provide one of the key experi- 
mental constraints on the theory of ultra-dense matter (e.g. Refs. [U [2]). 
The masses measured in the pulsar binaries are clustered around the value 
1.4 M & and have been consider as "canonical" for a long time. However, 
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in recent years mounting evidence emerged in favor of substantially heav- 
ier neutron stars with M < 2Mq. In particular, the recent discovery of a 
compact star with a mass of 1.97 M & measured through the Shapiro delay 
provides an observationally "clean" lower bound on the maximum mass of 
a compact star [3]. 

On the theoretical side it is now well-established that the emergence of 
new degrees of freedom at high densities softens the EOS of matter. For 
example, allowing for hyperons can reduce the maximum mass of a sequence 
of compact stars below the canonical mass of 1.4 M . A similar reduction 
may occur if a deconfinement to quark matter takes place, although the 
softening of the EOS in this case is less dramatic. Thus, the observation of 
2Mq mass neutron star is evidence that the ultra-dense matter in neutron 
stars cannot be soft, i.e., agents that will substantially soften the EOS 
are potentially excluded. Thus, one of the outstanding challenges in the 
theory of compact objects is to exploit fully the consequences of this recent 
observation and, for example, to provide new EOSs that are capable to 
produce compact objects as massive as the millisecond pulsar J1614-2230. 

The cooling of neutron stars provides another channel on the properties 
of dense matter, which is sensitive to the composition of matter and weak 
interactions therein. The cooling of compact stars can be divided roughly in 
the following phases. After the initial non-isothermal phase of rapid cooling 
from temperatures T ~ 50 MeV down to 0.1 MeV, a neutron star settles in 
a thermal quasi-equilibrium state which evolves slowly over the time scales 
10 3 - 10 5 yr down to temperatures T ~ 0.01 MeV [HE]. In this latter 
phase the core of the star is isothermal and the temperature gradients are 
concentrated in the envelope. The cooling rate of the star during this pe- 
riod is determined by the processes of neutrino emission from dense matter, 
whereby the neutrinos, once produced, leave the star without further inter- 
actions. The understanding of the cooling processes that take place during 
this neutrino radiation era is crucial for the interpretation of the data on 
surface temperatures of neutron stars. While the long term features of 
the thermal evolution of neutron stars are insensitive to the non-isothermal 
cooling stage, the subsequent route in the temperature versus time plane 
strongly depends on the emissivity of matter during the neutrino-cooling 
era. Thus, another outstanding question raised by the work of Ref. [3] is 
how the massive neutron stars, featuring quark matter, cool? The recent 
observation of the substantial change in the temperature of the neutron 
star in Cas A poses a further challenge for the theory to explain drastic 
short-term drop in the temperature of this neutron star [6l [7] . Presently, 
consistent calculations of the cooling of massive compact stars are virtu- 
ally absent. First steps in this direction have been taken by a number of 
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A robust feature of cold quark matter is its color superconductivity \\.2\ 
113] , Unfortunately, consistent, realistic simulations of cooling of compact 
stars featuring quark matter are far from trivial because of the complex 
phase structure of the quark matter at low temperatures and the substantial 
effort that is needed to understand transport and weak interactions in this 
matter. The study of neutrino emissivities and thermal conductivity of 
quark matter in the color superconducting state is in its beginning. In 
some cases, even crude estimates are not available for these quantities and 
a pressing requirements for any realistic simulation of thermal evolution of 
massive compact stars is the development of the knowledge of transport 
coefficients of these phases. 

2. Two solar-mass compact stars with hyperons and color 
superconducting quarks 

The existence of hybrid stars with two solar masses was predicted in a 
number of models, including those based on the MIT model and NJL models 
of quark matter and its superconductivity (for reviews see Refs. [121 03] )• 
Our recent study based on relativistic hypernuclear Lagrangians predicts 
stiff hypernuclear EOS above saturation density [T3]. This enables one to 
construct stable configurations with masses equal and above the measured 
1.97 solar-mass star. The resulting configuration have "exotic" matter in 
their interiors in the form of hyperons and quark matter, of which the quark 
matter is color-superconducting in the two-flavor 2SC and/or three-flavor 
CFL phases. 

Here we briefly describe the set-up and main results of Ref. [14]. The 
nuclear EOS, as is well known, can be constructed starting from a number 
of principles, see, e.g., [UI2]. In Ref. |14] a number of relativistic mean- field 
models were employed to model the low density nuclear matter. As is well- 
known, these models are fitted to the bulk properties of nuclear matter and 
hypernuclear data to describe the baryonic octet and its interactions |15t ll6|. 
The underlying Lagrangian is given by 

£>b = ^2 ^Bh^iidfj, - g^B^fi - l^g P BT ■ pp) - (m B - g a B(y)]^B 

B 

- h>m N {g aN af - ^c{g aN a) A + ^ ^a(*7 M ^ - m x )^x - -F^F^, 

(1) 

where the S-sum is over the baryonic octet B = p, n, A, S 1 * 1 ' , ips are 
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the corresponding Dirac fields, whose interactions are mediated by the a 
scalar, isoscalar-vector and isovector-vector meson fields. The next-to- 
last term in Eq. (JTJ) is the Dirac Lagrangian of leptons, Fn U is the energy and 
momentum tensor of the electromagnetic field. The parameters in Eq. (pQ) 
correspond to the NL3 parametrization [T7J. Computations were made also 
with the GM3 parameterization [15], with the result that hyperonic matter 
cannot be accommodated within this model. The choice of this specific 
parametrization was made because the nucleonic matter has the stiffest 
EOS compatible with the nuclear phenomenology. The mean-field pressure 
of the (hyper) nuclear matter can be obtained from Eq. ([1]) in the standard 
fashion [1]. 

The high-density quark matter was described in Ref. [2] by an NJL 
Lagrangian, which is extended to include the t' Hooft interaction term (oc 
K) and the vector interactions (oc Gy) [18] 

8 

C Q = VK^m - rh)^ + GV(V^ 7 V) 2 + G s ^[(V^) 2 + (^X a ^) 2 } 

a=0 

+ G D ^[C^75 e a/37 ea bc (^)^][(^c)^75 e p<J7 6^] 

- K {det/$(l + + det/$(l - 75)^]} , (2) 

where the quark spinor fields ip^ carry color a = r,g,b and flavor (a = u, d, s) 
indices, the matrix of quark current masses is given by fh = diagj(m M , m<i, m s 
\ a with a = 1,...,8 are the Gell-Mann matrices in the color space, and 
Ao = (2/3) If . Here Gs and Go are the couplings in the scalar and di-quark 
channels. The charge conjugate spinors are defined as ipc = Cip T and 
rpc = ^ T C, where C = i 7 2 7° is the charge conjugation matrix. The parti- 
tion function of the system can be evaluated for the Lagrangian ([2]) neglect- 
ing the fluctuations beyond the mean-field |18j . To do so, one linearizes the 
interaction term keeping the di-quark correlations A c oc (^c)^75 e ^abcV's 
and quark-anti-quark correlations a a oc ipa^Pa- The pressure derived from 
(El) reads 



1 18 f A ( 
P = T^Ylj ^ 2 jhl + 2Tln 



1 + exp 



N 
T 



AKa u a d <?s 



3 3 



- i^Ei A c i 2 - 2 ^E cj ' + i^( 2w o 2 + ^)+ E pi-po-b*, 

c=l a=l l=e-,fj,- 

(3) 

where T is the temperature, are the quasiparticle spectra of quarks, uiq = 
G v {QM\iptip u + ^ d i/j d \QM) and cfo = 2Gy{QM\i/jtil> s \QM) are the mean 
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Fig. 1. Left panel. Mass vs radius for configurations with quark-hadron transition 
density p tr — 3.4p for four values of vector coupling Gy/Gs — 0, 0.2, 0.4, 0.6. The 
purely hadronic sequence (i.e. the sequence that includes nucleons and hyperons) 
is shown by black solid line. The dashed lines and the gray solid lines show the 
branches where the 2SC and CFL quark phases are present. The filled circles mark 
the maximum masses of the sequences. The horizontal line shows the largest mass 
measurement to date [3] . Right panel. Properties of the stars as a function the free 
parameters Gy and ptr- The solid lines show the maximum mass configurations 
realized for the pair of parameters Gy and ptr- The dashed curves show the amount 
of CFL matter in the configurations via the ratio S = Rcfl/R, where Rcfl is the 
radius of the CFL core, R is the star radius. The parameter space to the right from 
(5 = line produces CFL stars. The parameter space below the dashed-dotted 5.1 
curve corresponds to stars containing 2SC matter. 



field expectation values of the vector mesons oj and 4> in quark matter, pi 
is lepton pressure, po is the vacuum pressure and B* is an effective bag 
constant. The quark chemical potentials are modified by the vector fields 
as follows fb* = diagf(fi u —uJo,fid — ^o,IJ's — <fto)- The numerical values of the 
parameters of the Lagrangians ([1]) and ([2]) are quoted in Ref. [14] . The mass- 
radius relationship for massive stars constructed on the basis of the EOSs 
described above is shown in the left panel of Fig. [1] together with the largest 
mass measurement to date M = 1.97 ±0.04/14© [3]. Masses above the lower 
bound on the maximum mass are obtained for purely hadronic stars; this 
feature is prerequisite for finding similar stars with quark phases. Evidently 
only for high values of vector coupling Gy one finds stable stars that contain 
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(at the bifurcation from the hadronic sequence) the 2SC phase, which are 
followed by stars that additionally contain the CFL phase (for higher central 
densities). Thus, we find that the stable branch of the sequence contains 
stars with quark matter in the 2SC and CFL phases. 

The right panel of Fig. Q] shows the changes in the masses and composi- 
tion of compact stars as the parameters of the model Gy and are varied. 
First, it shows the tracks of constant maximum mass compact stars within 
the parameter space. The decrease of maximum masses with increasing vec- 
tor coupling reflects the fact that non-zero vector coupling stiffens the EOS. 
In other words, to obtain a given maximum mass one can admit a small 
amount of soft quark matter with vanishing vector coupling by choosing a 
high transition density; the same result is obtained with a low transition 
density, but strong vector coupling, i.e., a stiffer quark EOS. For low tran- 
sition densities one finds 2SC matter in stars, which means that weaker 
vector couplings slightly disfavor 2SC matter. Substantial CFL cores ap- 
pear in configurations for strong vector coupling and almost independent 
of the transition density (nearly vertical dashed lines with 5 ~ 0.1 in Fig. 
[H right panel). Note that for a high transition density there is a direct 
transition from hyper-nuclear to the CFL phase. For transition densities 
blow 3.5/>o a 2SC layer emerges that separates these phases. On the other 
hand, weak vector couplings and low transition densities produce stars with 
a 2SC phase only. 



3. Thermal evolution of massive stars 

New features arise in the cooling behaviour of compact stars with the 
onset of quark matter in sufficiently high-mass models. The recent de- 
velopment of sequences of stable massive hybrid stars with realistic input 
EOS \19\ \20 \ [21] allows us to model the thermal evolution of compact stars 
containing quark cores. These sequences of stable stars permit a transition 
from hadronic to quark matter in massive stars (M > I.85M0) with the 
maximal mass of the sequence ~ 2M . In Ref. [11] quark matter of light u 
and d quarks was assumed in beta equilibrium with electrons. The pairing 
among the u and d quarks occurs in two channels: the red-green quarks 
are paired in a condensate with gaps of the order of the electron chemi- 
cal potential; the blue quarks are paired with (smaller) gaps of order of 
keV, which is comparable to core temperature during the neutrino-cooling 
epoch [221 1231 [2^] . For the red-green condensate, a parameterization of neu- 
trino emissivity was chosen in terms of the gaplessness parameter £ = A/<5/x, 
where A is the pairing gap in the red- green channel, 5p is the shift in the 
chemical potentials of the u and d quarks The magnitude of the gap 
in the spectrum of blue quarks was treated as a free parameter. 
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The stellar models were evolved in time, with the input described above, 
to obtain the temperature evolution of the isothermal interior. The interior 
of a star becomes isothermal for timescales t > 100 yr, which are required 
to dissolve temperature gradients by thermal conduction. Unless the ini- 
tial temperature of the core is chosen too low, the cooling tracks exit the 
non-isothermal phase and settle at a temperature predicted by the balance 
of the dominant neutrino emission and the specific heat of the core at the 
exit temperature. The low-density envelope maintains substantial temper- 
ature gradients throughout the entire evolution; the temperature drops by 
about 2 orders of magnitude within this envelope. The isothermal-interior 
approximation relies further on the fact that the details of the tempera- 
ture gradients within the envelope are unimportant if we are interested only 
in the surface temperature T s . Models of the envelopes predict the scal- 
ing T 4 = g s h(T), where g s is the surface gravity, and h is some function 
which depends on T, the opacity of crustal material, and its EOS. The fitted 
formula T 8 = 1.288(T 4 6 /c/ s i4) a455 [26j is commonly used. 

In the isothermal-interior approximation, the parabolic differential equa- 
tion for the temperature reduces to an ordinary differential equation, 



where L v and L 7 are the neutrino and photon luminosities, Cy is the spe- 
cific heat of the core, and the heating processes, which could be important 
in the photon cooling era, are neglected, i.e., H(T) = (see Refs. [27] 128] 
for a summary of these processes and their effect on the evolution). The 
results of integration of Eq. Q are shown in Fig. [51 where we display the 
dependence of the (redshifted) surface temperature on time. Each panel of 
Fig. [2] contains cooling tracks for the same set of four models with central 
densities 5.1, 10.8, 11.8, 21.0 in units of 10 14 g cm~ 3 . The cooling tracks for 
the purely hadronic model (solid lines) are the same in all four panels. The 
panels differ in the values of micro-physics parameters, which characterize 
the pairing pattern in quark matter. Specifically, the two panels in the left 
column correspond to the case where the blue-quark pairing is negligible 
(i.e., the pairing is on a scale much smaller than the smallest energy scale 
involved, typically the core temperature). The two panels in the right col- 
umn correspond to the case where the gap for blue quarks is large, = 0.1 
MeV. The panels in the upper and lower rows are distinguished by the value 
of the C parameter. [We use the values £ = 0.9 (upper row) and £ = 1.1 
(lower row)]. It can be seen that (i) the neutrino-cooling is slow for hadronic 
stars and becomes increasingly fast with an increase of the size of the quark 
core, in those scenarios where there are unpaired quarks or gapless excita- 
tions in the superconducting quark phase. The temperature scatter of the 




L„(T)-L^(T S ) + H(T) 



(4) 
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Fig. 2. Time evolution of the surface temperature of four models with central 
densities 5.1 (solid line), 10.8 (long-dashed line), 11.8 (short-dashed line), 21.0 
(dotted line) given in units of 10 14 g cm -3 . For the observational data sec Ref. [llj . 
The upper two panels correspond to cooling when the red-green condensate has 
C = 0.9, i.e., is not fully gapped; the lower panels correspond to ( = 1.1, i.e., the 
red-green condensate is fully gapped. The left two panels correspond to evolution 
with negligible blue-quark pairing (A& = 0); the right two panels show the evolution 
for large blue pairing A& = 0.1 MeV. 



cooling curves in the neutrino cooling era is significant and can explain the 
observed variations in the surface temperature data of same age neutron 
stars, (ii) If quarks of all colors have gapped Fermi surfaces, the neutrino 
cooling shuts off early, below the pairing temperature of blue quarks; in 
this case, the temperature spread of the cooling curves is not as significant 
as in the fast cooling scenarios, (hi) As the stars evolve into the photon 
cooling stage the temperature distribution is inverted, i.e., those stars that 
were cooler in the neutrino-cooling era are hotter during the photon cooling 
stage. 
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4. Perspectives 

The physics of massive compact stars poses a number of interrelated 
questions/challenges. The first issue is the EOS of matter, including the 
quark degrees of freedom and their color superconductivity. The strangeness 
degrees of freedom including hypernuclear matter and three-flavor quark 
matter need to be further explored building, e.g., upon the work of Ref. [14] . 
The equilibrium and stability of massive compact objects, constructed from 
these EOSs, should be studied including rapid rotations and oscillations. 
Secondly, we need a better understanding of the weak interaction rates in 
quark and (hyper) nuclear matter, which are required input in cooling sim- 
ulations of compact stars. Thirdly, the transport coefficients of dense color 
superconducting quark matter, such as the thermal conductivity, are needed 
for modelling an array of phenomena, which include thermal evolution, mag- 
netic evolution, r-modes etc. 
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